Spin waves in the spiral phase of a doped antiferromagnet: a strong-coupling approach 
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We study spin fluctuations in the spiral phase of the two-dimensional Hubbard model at low dop- 
ing on the basis of the spin-particle-hole coherent-state path integral. In the strong correlation limit, 
we obtain an analytical expression of the spin-wave excitations over the entire Brillouin zone except 
in the vicinity of q = 0. We discuss the validity of the Hartree-Fock and random-phase approxima- 
tions in the strong-coupling limit, and compare our results with previous numerical and analytical 
calculations. Although the spiral phase is unstable, as shown by a negative mean-field compressibil- 
ity and the presence of imaginary spin-fluctuation modes, we expect the short-wavelength fluctuation 
modes (with real energies) to survive in the actual ground-state of the system. 
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I. INTRODUCTION 



Despite a lot of theoretical efforts, there is still no sat- 
isfying description of the ground-state and the low-lying 
excitations of a doped two-dimensional (2D) antiferro- 
magnet. One of the simplest (realistic) models describ- 
ing such a system is the Hubbard model. In the strong 
correlation limit, it reduces to the Heisenberg model at 
half-filling. For a 2D square lattice, the ground-state is 
known to be antiferromagnetic, and the low-lying excita- 
tions (spin-waves) are well understood. However, away 
from half-filling, the coexistence of local moments and 
itinerant charge carriers raises difficulties that have not 
been overcome so far. 

In the framework of the Heisenberg model, perturba- 
tion theory around a broken-spin-symmetry ground-state 
(from linear spin-wave analysis to renormalization-group 
approach) has proven to be very successful. Formally, 
this semiclassical approach corresponds to a 1/S expan- 
sion, where S is the size of the localized spins. It is very 
natural to follow the same line of approach in itinerant 
spin-i fermion systems by promoting the localized mo- 
ments from spin-i to spin-S (with S 3> 1). The spin de- 
grees of freedom are then treated semiclassically, whereas 
the itinerant fermions are considered quantum mechani- 
cally. 

The semiclassical approach is realized for instance 
in the Hartree-Fock or random-phase (RPA) approxi- 
mations. For the Hubbard model, the (homogeneous) 
HartreepEock theoryEriiJ and the slave-boson mean-field 
theoryaO predict a spiral magnetic order at strong cou- 
pling in agreement withjaaean-field (semiclassical) theo- 
ries of the J modell3il3 as first shown by Schraiman 
and Siggia.EZI This (homogeneous) spiral phase turns out 
to be unstable, as shown by a negative compressibil- 
ity. The instability also manifeata-.itself by the presence 
of imaginary spin-wave modes B'BHii3 Short-range spiral 
order,EZI, phase separation, coexisting spitt-H-and charge- 
density waves (domain- wall fesriaation) ,EJ'BEra and for- 
mation of local spin polaronsEjEil have been proposed 



as a possible alternative to the spiral phase. Although 
phase separation is likely to be suppressed by long-range 
Coulomb interaction, it is not clear whether the latt 
can stabilize the homogeneous spiral magnetic order! 
In any case, we expect the short-range spiral order to 
survive in the actual ground-state of the system. The 
study of (short-wavelength) spin-wave modes around the 
mean-field spiral order can then be seen as a way to ac- 
cess this short-range order, even though the mean-field 
theory incorrectly predicts a long-range spiral order. On 
the experimental side, strong antiferromagnetic fluctua- 
tions have been observed in the normal phase of high-T c 
copper oxides. In particular, inelastic neutron scattering 
data have revealed noticeable incommensurate fluctua- 
tions at finite doping E3 

In this paper, we reconsider the semiclassical limit of 
the 2D Hubbard model on the basis-pf the spin-particle- 
hole coherent-state path integral.caca The latter has been 
designed to study strongly correlated fermion systems 
where itinerant charge carriers interact with local mo- 
ments. As shown below, it provides a very natural frame- 
work to derive the spin-wave modes around the mean- 
field spiral order. 

The outline of the paper is as follows. In Sec. H we in- 
troduce the effective action of the Hubbard model based 
on the spin-particle-hole coherent-state path integral. Af- 
ter performing a saddle-point approximation on the spin 
variables corresponding to a (homogeneous) spiral order, 
we calculate the hole Green's function and the free energy 
within a t/U expansion. We discuss the reason why the 
Hartree-Fock theory turns out to be correct in the strong 
correlation limit at low doping. In Sec. Ill, we compute 
the spin-wave excitations. We point out that a correct 
description of the three Goldstone modes of the spiral 
phase (at q = 0,±Q, where Q is the wave-vector of the 
spin modulation) follows from a consistent calculation of 
vertex corrections and self-energy terms in the fermion 
Green's functions. Finally, we obtain an analytical ex- 
pression of the spin-wave excitations valid over the entire 
Brillouin zone, except in the vicinity of q = 0. Wer 
pare our analytical result with previous numericalE 



1 



and analytical calculations 
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II. MEAN-FIELD THEORY 

We consider a bipartite 2D lattice with N sites. The 
Hubbard model is defined by 
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' a + h.c.) + U n r i-n r ;, 



(2.1) 



where c rCT is a fermionic operator for a er-spin particle at 
site r (cr =t,|)) fi ra = cJg-Cro-, and (r, r') denotes nearest 
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neighbors. We denote by fi the chemical potential, j3 the 
inverse temperature, and n — 1 — x the mean number of 
particles per site. We consider only hole doping (x > 0) 
and set Ti = ks = 1 throughout the paper. All results 
are obtained in the zero-temperature limit (T — * 0). 

The spin-particle-hole coherent-state path integral 
formulation of the Hubbard model was derived in 
Refs. |23| , p4] . This approach is based on the introduction 
of spin-particle-hole coherent states which generalize the 
spin-i coherent states by allowing the creation of a hole 
or an additional particle. In the strong-coupling limit 
U 3> t, the effective action 5[7*,7;f2] of the Hubbard 
model is given by 
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dr 7r*T ( 9 r ~ V + ^r)7rT - E 

r,r' 

dT^idr-fl + U-A^)^ 
dr (7rT^T,r'l7r'l + C.C.J 



^ T 7r*T*rT,r'T7r'T 

dr 7rlil,r'l7r'l 
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+ E / dn dT 2 dT 3 dT4r\\ ]Ti (n,T2;r 3! r4)7; T (Ti)7* i (r 2 )7 r |(T4)7 rT (T3). 



(2.2) 



$7 is a unit vector field which gives the direction of 
the local moments at the singly-occupied sites. 7j and 
7j_ are Grassmann variables describing particles propa- 
gating in the lower (LHB) and upper (UHB) Hubbard 
bands, respectively. A® — (tt r \Cl r ) is a Berry phase 
term (|O r ) = 9 T |f2 r )). The intersite hopping matrix 
t rr > = Rlt rr iR r > depends on f2 via the SU(2)/U(1) matrix 



R, 



(2.3) 



which rotates the unit vector z to f2 r = fl(6 r , <Pr)- Here 
8 r and ip r are the polar angles determining the direc- 
tion of Jl r . The choice of ip r is free and corresponds to 
a "gauge" freedom. In the Hubbard model, t rr / = t if 
r and r' are nearest neighbors and vanishes otherwise. 
The quartic term of the action ( [2^ ) is determined by the 
two-particle atomic vertex r . 

In this section, we make a saddle-point approximation 
on the spin variables f2 r . As discussed in Ref. ^4], this 
approach can be justified by taking a "large-5"' semi- 
classical limit, which consists in promoting the spin-^ 
coherent states describing singly occupied sites to spin-S 1 



coherent states. In the limit S — * oo, the Berry phase 
term suppresses quantum fluctuations of f2. The spin 
variables become classical and do not fluctuate at zero 
temperature. 

A broken-symmetry ground-state corresponding to a 
(coplanar) spiral order is defined by the "classical" con- 
figuration 
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<pf = Q r. 



(2.4) 



We consider only the diagonal spiral phase (Q = (Q,Q)), 
which is known, to be the most stable one in the strong- 
coupling limitfl Note that the antiferromagnetic (Q = 
(tt, 7r)) and ferromagnetic (Q = 0) phase s a re special 
cases of the spiral order defined by Eq. (2.4). In the 
gauge ip r = 0, the hopping matrix i™, depends only on 
the difference r — r': 



t cl , = t 



-§Q-(r-r')<r x 



(2.5) 



This yields the following saddle-point action: 



U5 a i + 2t cos y cos/c l/ )7k (T («w) + ^2 7LM 2t sin y sinfc 1 , 7kff (^) 
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(2.6) 
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where 7k CT is the Fourier transform of 7 r(T , and k runs 
over the entire Brillouin zone, T 1 ^ ^(iwi, iu>2i ^3(^4)) = 
iu>2 — i^>3 — U, where W4 = u>i+uj2~ll>3 is fixed by energy 



conservation. In Eq. (2J3) and below, there is an implicit 
sum over v = x, y. Without the quartic term, the ac- 
tion (2.i) is similar to that obtained in Refs. |||| within 
a large-t/ Hartree-Fock approximation where the SU(2) 
spin-rotation invariance is maintained by introducing a 
fluctuating spin-quantization axis in the functional inte- 
gral. We show below that the quartic term can indeed 
be neglected in the low- doping limit. 

For later convenience, we introduce the parameter 



cos — > 



(2.7) 



which determines the spiral pitch. By changing Q in 
Q + 2ir, it is always possible to choose the sign of p posi- 
tive, p varies from in the antiferromagnetic phase to 1 
in the ferromagnetic phase. 
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FIG. 1. Definition of the various symbols appearing in the 
Feynman diagrams. 
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FIG. 2. Diagrammatic representation of the LHB 
self-energy at order 0(t/U). Slashed dashed lines indicate 
interband transitions (i%g-). 



We are now in a position to apply the strong-coupling 
perturbation theory described in Ref. |24|. We first con- 
sider the LHB Green's function. To zeroth order in t/U, 
i C Q S is ignored, since it corresponds to interband transi- 
tions. As shown in Ref. ^4|, T n does not affect the LHB 
Green's function to this order. We thus obtain 



£/j (k, iuj) = (ioj + fx — 2tp cos k v ) 1 



(2.8) 



as the propagator of the 7j field. Corrections to are 
taken into account by introducing a self-energy E, i.e. 
Q^ 1 = 1 — E. To order 0(t/U), there are three con- 
tributions to E shown in Fig. §B [The symbols used in 
the Feynman diagrams are defined in Fig. g.] The last 
contribution (Fig. ||c) vanishes due to the sum over the 
internal momentum in the loop. We show at the end of 
this section that the second contribution (Fig. |^b) can be 
neglected in the low-doping limit. The first contribution 
(Fig. §a) gives 



Z(k,iu) = 4< 2 (l-p 2 )(sinfc I ,) 2 G{ 0) M 
~ -J(l -p 2 )(smk u ) 2 , 



(2.9) 



where J = 4t 2 /U. The last line of Eq. (gj) is ob- 
tained by approximating the UHB atomic Green's func- 
tion Gf\iLu) = (ioj + n- U)- 1 by -l/U. We therefore 
obtain the following expression for the propagator of the 
7 field: 



Qf(k,iu}) = (iuj - ek) , 

2tpcosk ly - J(l - p 2 )(smk v ) 2 . 



(2.10) 
(2.11) 



It should be noted that Q-t is not the LHB Green's func- 



tion of the original fermions (i.e. the c field in Eq. (2.1)). 
The latter is given by (see Eq. (3.32) in Ref. ^4|) 

gJ HB (rr,rV) = (i? r (r)) CTT ^(rr, r'r')(R' r (r')) U 

--jcrQ-(r-r') 

= o ST(rr,rV). (2.12) 



In Fourier space, this gives 



1 



^(k^) = -£ T (k+-Q, iW ). 



(2.13) 



The quasi-particle pole has a residue equal to 1/2, as ex- 
pected since <7^ HB is the Green's function projected onto 
the LHB. Moreover, the dispersion law of the original 
fermions is shifted by ±Q/2 with respect to that of the 
7f parti cles. 



Eqs. (2.11) and (2.13) agrees with the large-C/ expan- 
sion of the Hartree-Fock result. There are two contri- 
butions to the energy of the particles in the LHB. 
The first is due to inter-sublattice hopping processes and 
gives a band of width 8tp. The second comes from 
intra-sublattice hopping processes which occur via vir- 
tual transitions to the UHB. The resulting bandwidth 
is of order J(l — p 2 ). Anticipating that p ~ (U/2t)x 
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[Eq. (2.21)], inter- (intra)-sublattice hopping processes we obtain a parabolic band with anisotropic effective 



dominate when J(l — p ) <C tp, i.e t /U <C a; (ip <C masses 



J(l -_p 2 ), i.e. £ < t 2 /U 2 ). From Eq. ( [2.1l| ), we readily 
see the instability of the antiferromagnetic phase against 
the formation of a spiral phase in the presence of holes. 
When p is finite, the kinetic energy gain due to inter- 
sublattice hopping (or order p) dominates over the en- 
ergy loss due to intra-sublattice hopping (of order p 2 ), 
thus stabilizing a non-linear spin order. [Note that the 
energy loss due to intra-sublattice hopping can also be 
seen as a loss of exchange energy] 




-1 




-3-2-10 1 2 3 

FIG. 3. Dispersion law e k in the LHB for t 2 /U 2 < x 
(U/t = 10 and x = 0.15). 

Since p > 0, the top of the LHB (for the 7 particles) 
is located in the middle of the Brillouin zon e (k = 0). In 
the vicinity of k = 0, the dispersion law (2.11) can be 
approximated as 

e k ~ -/i + Atp - tpk 2 - J(l - p 2 )(k x + k y ) 2 . (2.14) 

Introducing the momentum K defined by 

—k x -\- k v 



_ k X + ky 

X V2 ' V ~ V2 



(2.15) 



2m 7 



Kr, 



2m 



u 



— = 2tp + 4J(l -p 2 ), — =2tp. 
m x m y 



(2.16) 



The dispersion law ek is sh own in Figs. ^ and (we have 
used p ~ (U/2t)x, see Eq. (^.2l|)). When t 2 /lP < x, the 
anisotropy of the effective masses is weak and the Fermi 
surface is almost circular. The minima of the LHB lie 
at the corners of the Brillouin zone (Fig. |). With de- 
creasing x, the Fermi surface becomes elliptical, and the 
minima of the LHB move away from the corners of the 
Brillouin zone (Fig. 0) . 
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FIG. 4. Dispersion law e k in the LHB for x ~ t 2 /U 2 
(U/t = 10 and x = 0.04). 

For T = 0, the free energy (per site) F(^) is given by 
(see Appendix ^) 



-ek)ek 
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-ix- J(l-P 2 ) 



47T 



^m x m y . (2.17) 



The chemical potential is obtained from n = 1 — x = 
-dF/dn?3 i.e. 



ji = Up- 



2ffX 



(2.18) 



From Eqs. (2.17) and (2.18), we deduce the energy 
E{x) = F + iKT^ x): 



E = -ttpx - J(l - p 2 ) 



(2.19) 



The spiral pitch is obtained by minimizing the energy E 
at fixed hole doping x, i.e. dE/dp = 0. This gives the 
following equation for p: 



2.2 : v , ,„n _p2^i/2 

2„ I QJ- 7/1 o„2\ 



(4te - 2Jp)[4rp 2 + 8tJp(l -p 2 )] 



— [8ty + 8tJ(l - 3p 



In the limit x -C 1, we fin 



P 



27" 



(2.20) 



(2.21) 



At half-filling, the spin configuration is antiferromagnetic 
(p = 0). As soon as holes are introduced into th e sys - 
tem, a spiral phase is stabilized. Since p < 1, Eq. (2.21) 
holds only when x < x c — 2t/U. When x > x c , the 
minimum of the energy is reached for the ferromagnetic 
state (p = 1). 



From Eqs. ( 2.18 ) and ( 2.21 ), we deduce that the com- 
pressibility k = dfi/ dfi — —1/2U is negative in the spiral 
phase. This signals an instability towards phase sepa- 
ration. The presence of imaginary spin-wave modes at 
finite momenta (see Sec. [II) suggests that other types of 
ground-states could also be stabilized. Zhou and Schulz 
have proposed a charge- density wave (coexisting with 
some kind of magnetic order) .013 This possibility has also 
been studicd-in the framework of the Schraiman and Sig- 
gia's model.ll3 The ferromagnetic phase, which has a pos- 
itive compressibility k ~ l/47rf, is found to be stable. 
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FIG. 5. Lowest-order correction F' to the free energy F 
due to the two-particle vertex r . 

Let us now consider the self-energy shown in Fig. ^d. 
Although its calculation presents no difficulty, the deter- 
mination of its contribution to the free energy requires 
to perform a coupling constant integration in order to 
ensure a consistent description of the system thermody- 
namics. Instead of considering X, we evaluate the lowest 
order correction F' to the free energy F due to the two- 
particle vertex T n (Fig. ||). This is sufficient to show that 



r 11 can be ignored in the low-doping limit (i « 1). As 
shown in Appendix H one finds 

F' = xJ(l — p 2 ). (2.22) 

F' can be sim ply taken into account by replacing J(l— p 2 ) 



inEq. ( [2.17 ) by J(l— p 2 )(l — x). This modification can be 
ignored for x <C 1. We conclude that T n can be ignored 
in the low-doping limit. 

The fact that we can neglect the quartic term of the 
action £[7*, 7; ft] justifies the functional integral formu- 
lation of Refs-Eip. Ji is also the reason why a standard 
RPA analysislffl'ErEll gives correct results in the strong- 
coupling limit. Note however that this conclusion holds 
only in the low- doping limit [x <C 1). For an arbitrary 
doping, the quartic term of the action (2.6) cannot be 
neglected, and the RPA approach breaks down since it 
misses processes of order 0(t/U). It should also be noted 
that the UHB Green's function Qi obt aine d by diagonal- 
izing the quadratic part of the action (2J3) (which is the 
starting point of the RPA approach) is meaningless in 
the strong-coupling regime since it does not correspond 
to the leading term within a t/U expansion (see Ref. |24| 
for a detailed discussion of this point). 



III. SPIN- WAVE MODES 
A. Effective action 

In this section, we derive the effective action for the 
fluctuations 8Q r = fl r — of the spin variables around 
their saddle-point value. The calculation is performed 
within a perturbative expansion in t/U . j—. 

We parametrize the fluctuations 8fl r bycll 



p r = i<50 r • r , q T = Sft r ■ ip r , 



(3.1) 



where 8 r = z and Cp v = {— sin(Q-r), cos(Q-r), 0). q T (p r ) 
corresponds to fluctuations in (out of) the spiral plane. 
We easily obtain 



Pr 



cos 9 r 



fir 



'■Pr 



Q r, 



(3.2) 



where the last equality holds for small fluctuations. The 
Berry phase term can be expressed in terms of the vari- 
ables p, q: 
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-gCPrSr - Pa- 



using Eqs. (2.3) and (3.2), we deduce 
Si 1 = i , -i cl , =t 

where 



A rr , -b; 

£rr< A* 



(3.3) 



(3.4) 
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J 



A rr > = cos 



B rr > = COS 



Q (r r') 



Q (r r') 



Q (r r') 



Q (r r') 



{Pr-Pr') 2 , 1 



+ +Pr')(«r ~ <7r' ) 



<7r - <7r 



(<7r - g r ') 5 



Pr-Pv' - -j[q r 



% 1 Z 

^{Pv+Pr') 2 - ^{jPr ~ Pr>){qr - q r >) + g(<7r - <7r') 2 



(3.5) 



(3.6) 



to quadratic order in p, q. 

Taking into account spin fluctuations, the action can 
be written as 

s[ 7 *,7;tt] = W,7] ~E J *nl*i**T*' 

r,r ; 

+ E ff / drA° rl * ralra , (3.7) 



S B = - x) ^ / dr(p r q r - p r q r ) 



(3.13) 



5b is the standard expression for the Berry phase term 
of localized spins, with a reduction factor 1 — x due to 
doping. 



2. First-order cumulant Si 



where 7 r = (7 r f , 7ri) T - The action of the spin degrees of 
freedom is obtained by integrating out the fermions: 



SM = / x>[ 7 ] e -s 



-S[7*>7iO] 



We write the first-order cumulant as S% = S[ + S" 
where S{ (S[') is of order 0(1) (0(t/U)). The corre- 
sponding Feynman diagrams are shown in Fig. We 
(3.8) find 



For small fluctuations around the classical configuration 
fl , it is sufficient to determine S[p, q] to quadratic order 
in p, q. We find 



Si 



where 



S\p,q] = S B [p,q] + Sx\p,q] + S 2 \p,q], 
S B=^2<7 / dr^(7 r * CT 7 r(T ) c i, 

rr' 7r' / cl j 



Si 

s 2 



(3.9) 

(3.10) 
(3.11) 



J dTSi rrtr , r g r (r' -r,r = 0-) 

r,r' 

■^«5t kT ,kT^(k,^), (3.14) 



S" = -J2 J dTbi ra ^ s Q Sa {v' - r,r = 0-) 



r,r ,cr 



k,o;.(7 

where we have introduced 



r,r' 

l 

2 



^ / dTidr 2 (7 r t 1 (5£ rir ' i (ri)7 ri 

ri,ri,r 2 ,4 
X 7 r t 2 ^r 2 r' (T2)7r')cl- 



6t 



kcr-k'tr' ~~ 



1 E' 



:- ikr+ik ' r 'SL 



(3.15) 



(3.16) 



(3.12) 



The Green's function £/j(k, iw) is given by Eq. ( 2.10 ). To 
first order in t/U, 



(• • -) c i means that the average is taken with the saddle- 
point action 5 c i[ 7 *, 7], only the connected part being con- 
sidered. Sb is a Berry phase term, Si a first-order cu- 
mulant, and S2 a second-order cumulant. Note that A 

does not contribute to the second-order cumulant, since where = sinfc„. 
it is of second order in p, q. 



1. Berry phase term Sb 
Since (7* .7 ro -)cl = (1 - x)5 a ^, we obtain 



Q,rAk.U-) ■■ 2tsin|x k g T (k,iu,)Gf ) (zu;) 
2t . Q 



(3.17) 




dt tt !! w %t 

FIG. 6. Diagrammatic representation of S[ and 5*" . 
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Performing the sums over the Matsubara frequencies 
and using 

St U ■ Q Y 

v Jy * 



t 



■ ql 2 + J^sL ] ( y k _ y k+q ) 



i^cos — ^p q g_ q (y k+ q-y k _ q ), (3.18) 
q 

4t . Q 



^ <5t k(T , kff = — sin — [bq| 2 (^k + A k+q ) 



k q | 2 



(A k — A k+q ) 



(3.19) 



As will be shown below, generating the diagrams for 
S\p, q] from the self-energy S ensures a proper descrip- 
tion of the Goldstone modes. A mere calculation of the 
action S[p, q] to order 0(t/U) would yield a contradiction 
with Goldstone theorem. 

(a) 




<% M uSttf Sttf „ 




we finally obtain 



2tp 
k,<j 

2J(1-P 2 ) 

2 



^ = -^E^- £ 0(>k-ik +q )(i^ 



5f 



1 9g I 



^0(-e k )A k [b 9 ~| 2 (A k + A k+q ) 

k,g 



Ki<7 



-(A k — ^ k + q ) 



(3.20) 



where — cos fc„. Here we use the notation q — (q, iu> v ) 
with uj v a bosonic Matsubara frequency. Pq and are 
the Fourier transforms of p r and q r with respect to space 
and time. 



3. Second-order cumulant So 



We decompose 52 into two contributions, S 2 and S 2 , 
which are represented diagrammatically in Fig. ^. S 2 in- 
volves particle-hole excitations in the LHB and describes 
the dynamical interaction between holes and spin fluctu- 
ations. This interaction can also be seen as an RKKY- 
type interaction between localized spins mediated by the 
mobile holes. E£l The first diagram of Fig. 0a gives a contri- 
bution which is O(l) in t/U. The other three diagrams 
take into account "vertex corrections" to the external 
"vertices" St^ appearing in the O(l) contribution. These 
vertex corrections can be systematically generated from 
the self-energy diagram shown in Fig. |]a. Replacing 
by i^g. + 8i aS in this diagram, one obtains the vertex cor- 
rections shown in Fig. ga. Including the latter into the 
O(l) contribution to S' 2 , we obtain all the diagrams con- 
tributing to S' 2 (Fig. Ra). This procedure also generates 
a diagram of order Sp (Fig. ||b), which gives rise to S 2 
(Fig. [7p). It should be noted here that S' 2 contains dia- 
grams (but not all of them) of order 0(t 2 /U 2 ). We also 
note that the diagrams contributing to the first-order cu- 
mulant Si (Fig. ^|) can be obtained in the same way, i.e. 
by considering the vertex corrections shown in Fig. pk. 



(b) 



St 



u 




FIG. 7. Diagrammatic representation of S'2 (a) and S'2 (b). 



(a) St 



It 



St 



n 



(b) Si 




FIG. 8. (a) "Vertex corrections" generated from the 
self-energy E (Fig. |a) by replacing i% by i% + 6t aS . (b) 
Diagram 0(5t 2 ) generated by the same procedure. 



Let us first consider S 2 . We find 



dTl dTzSt^u^ a(Tl)Sir 2 ay 2 cr(T2) 



x [^, T Gf « - r 2 ,n - r 2 )£ T (r 2 - n,T 2 - ti) 
+S ad g^r[ - r 2 ,n - T 2 )Gf\r' 2 - n,r 2 - n )l. (3.21) 



Since spin fluctuations have a characteristic energy scale 
J -C U, we can replace 8t(r 2 ) by St(ri) in Eq. (3.21). 
This allows to obtain 



sa = 



-y 

2U ^ 



dTSik'ak&Stkvk'a 



k ,kV 

x[S atl 9(-e k ) + 5^6(-ev)}. 



(3.22) 



To linear order in p, q, 
^*kT,k'j. 



-j= cos — \Y V - F k )(^Pk-k' + 2 9k - k ' 



N 



7 



HQ. .( i 

o*k|,k'T = —rrr cos — (^k' - >k)(^Pk-k' - 2 9k - k ' 



>N 2 
so that we finally obtain 
Jp 2 



N 



5>(-e k )(y k+q -r k ) 2 (|p,- 



k,g 



(3 23) * s ^ ne polarization bubble of the fermions in the LHB. 8t R 
is a rcnormalized vertex which takes care of the vertex 
corrections discussed above. To linear order in p, q 

2t 



(3.24) 



The sum of S\ and S% can be written in the form 

4 -(Aq-A, 



^ +#2 =E [|W! 2 (^ + A q ) + i^-(/- q - A q ) 



(3.25) 



where 



^ = jv E ^- £ x) [- 2 ^(^k - *k +q ) 

k 

+2.7(1 - p 2 )X 2 - V(Yk +q - Yk) 2 
* 2J(l-p 2 ) 



AT 



2_j ^(~ e k) A" k X k+q . 



(3.26) 



Let us now consider S' 2 . It can be written as 



S ? = ^E^k+ q T,kT(^)^kT,k+ q T(-*^) n ( k ; ?)> 

(3.27) 



k,g 



where 



n(k; g) = ~ E ( k > H^T ( k + q, «w + iu v ) 



i(_ ek ) _ g(- £k+q ) 
£k - £k+ q + ioj v 



(3.28) 



or k'T,kT 



sin — 

N 2 



Pk'-k(^kA k ,k' +-^k'Ak',k) 

(3.29) 



"2 ( i , k'-k(ArkAk ) k' — Ak'Ak',k) 



with 

A k , k ' = 
We finally obtain 



l-|cos|(F k -y k 0. 



(3.30) 



ft (9) 



ift-(g) 



Pq 



where we have introduced 
1 



S±(q) 
S+-(q) 



2N 



En(k;?)A±(k;q), 



(3.31) 



(3.32) 



i-E n ( k ^1A + (k;q)A_(k;q), (3.33) 



2N 



and 



■ Q, 



A±(k;q) = 2isin|-(X k A k , k+q ±X k+q A k+qik ). (3.34) 



The total action S = Sb + Si + S2 is then written as 



J 



hq + Aq + S+(q) 
x) 



Pq 
Qq 



^-±5 + _(ff) |(/i q -A q + S_(g)) 

The spin-wave excitations are obtained from the equation 

[hq + A q + S+ (q, w)] [/i q - A q + S_ (q, oj)] = [(1 - a;)w + S+_ (q, w)] 2 , 



(3.35) 



(3.36) 



where S(q, oj) = S(q , ioj v — > a; + i0 + ) is the retarded 
part of 5(9). Eq. ( |3.36| ) was first obtained in Refs. ||^. 



B. Goldstone modes 

In this section, we sho w th at the equation for the 
spin- wave excitations [Eq. ( 3.36 )] accounts for Goldstone 
modes at q = and q = ±Q in agreement with— the 
Goldstone theorem for a spiral spin arrangement Ball3 

A_(k;q), SL(q, oj) and S^ (q, oj) vanish at q = 0. 



Since /i q= o — A q= o, we conclude that oj = is solu- 
tion for q = 0. This mode satisfies p r = and therefore 
corresponds to a global rotation in the spiral plane. 
When q = Q, we have the properties 

h Q + A Q = i E 6 , (-£k)(ek+Q - ek), 

k 

S+(Q, oj = 0) = - i E ^(-ek)(ek+Q - e k ). (3.37) 



The last equality follows from A + (k;Q) = £k+Q — £k- 
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We have also used the fact that 5+(Q,oj = 0) is purely 
real since particle-hole excitations are gapped for q = Q 
at low doping (see Fig. |l^). Using A±(— k — q; q) = 

=FA±(k;q), one can also show that S-\ (q,w = 0) = 0. 

This result, together with /iq + Aq + 5+(Q, lo = 0) = 0, 
ensures the existence of a gapless mode for q = Q. A 
similar reasoning holds for q = — Q. Goldstone modes at 
±Q satisfy q r — and correspond to fluctuations out of 
the spiral plane. 



C. Analytical expression of the spin- wave excitations 

The long- wavelength snipp wave modes have been stud- 
ied in detail before. t9&Btffl They are observed only for 
very low doping. Above a certain hole concentration, 
these modes completely dissolve into the particle-hole ex- 
citation spectrum. From now on, we focus on the region 
q ^S> Kf(o), where Kf{q) — 0(y/x) is the "Fermi wave 
vector" in the direction specified by the angle a (see Ap- 
pendix [A|). In this momentum range, it is possible to 
derive an analytical expression of the spin-wave excita- 
tions. 

As shown below, the spiral phase sustains both real 
and imaginary modes. The latter signal an instability of 
the spiral phase. However, the (short-wavelength) modes 
with a real excitation energy are expected to survive in 
the actual ground-state of the system. Since they lie be- 
low the particle-hole excitation continuum (see Fig. 10), 

the imaginary part of S±(q, a;) and S-\ (q, lu) does not 

need to be considered. Furthermore, u) can be neglected 

against £k+q — £k in the calculation of S± and Since 

5-| (q, u = 0) = 0, the equation for the spin- wave modes 



J 



J 



2{\- P 2 ) + V(2-y q ) 2 + (i- P 2 ) 



then reduces to 



(1-x) 2 luI = [/i q + A q + S+(q)P c 



A q + S_(q)], 

(3.38) 



where S±(q) = S±(q, oj = 0). The summations over k 
involved in hq, A q and S±(q) are easily carried out in 
the low-doping limit. For an arbitrary function /(k), we 
have 



N 



k 



N 



where we have used (l/N) J3 k 0(ek) = x. In general, 
the sum over the entire Brillouin zone is easily evaluated. 
Eq. (3.3S) gives the leading contribution for x <C 1. Thus 
we find 



/i q = 2J(l-p 2 )-.V(2-Yq) 

+x[2tp(2-Yq) + Jp 2 (2- Yq) 
= 2J(1 -p 2 ) + xJp 2 (2~ Yq) 2 , 
A q = J(l-p 2 )F q , 

2tp 



S(q) = xU 2 (l-p 2 )- 



tq=0 



u 



(3.40) 
(3.41) 
(3.42) 

(Yq-2)! 2 , (3.43) 



where 5(q) = S+(q) = 5_(q). Eq. (|3.4]j ) has been ob- 
tained using p = (U /2t)x and is therefore no t val id in 
the ferromagnetic phase (x > x c ). From Eq. (3.38), we 
deduce 



*3 



(i - x(y q - 2)) 



21 2 



Y a 



,1-p 2 



Y a -2 -I 



1/2 



-[(i-p 2 W 



(3.44) 



At half- filling, Eq. ( |3.44 ) reproduces the spin- wave ex- 
citations of the Heisenberg antiferromagnet 



23 



cos q x + cos q y 



1/2 



(3.45) 



Thus Eq. ( p.44| ) turns out to be correct over the entire 
Brillouin zone when x — 0. However, away from half- 
filling, it is not valid in the vicinity of q = since it 
does not describe the long-wavelength Goldstone mode. 
Using (1 — p 2 )X±Q — p 2 {Y±q — 2) 2 , one easily verifies 
that u>±q 



0. Thus the analytical re sult (3.44) satis- 
fies the Goldstone theorem. From Eq. (3.44), we deduce 



that uiq is either real or purely imaginary. Figs. |9j and 
|l0| shows the real and imaginary parts of w q . We ob- 



tain a ver 
Refs. | 



rood agreement with the numerical results of 



r 



The existence of complex frequencies o; q in some re- 
gions of the Brillouin zone signals an instability of the 
spiral phase. The origin of this instability can be traced 
back by considering the static correlation functions (with 
the analytic continuation iio v — ► u>) 



(PqP*q)u=0 

(<Zg<z|)w=o 



1 



2{h c 



•5(q))' 



ft, - A q + 5(q) ' 



(3.46) 



where hq, A q and 5(q) are defined by Eqs. ( |3. 41 - 3. 43 ) 



While (PqP~} u=0 is always positive, (g r g?|) w=0 turns out 
to be negative where the excitation energy u Q be comes 
i magi nary (see Fig. |l^). [Compare Eqs. ( 3.38| ) and 
(3.46).] This implies a negative spin stiffness (and there- 
fore aHrnkistability) for spin fluctuations in the spiral 
planemB As shown in Refs. |,|, the imaginary modes 
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also lead to a negative charge susceptibility. This sug- 
gests that the actual ground-state could exhibit a charge- 
density wave (coexisting with spin order). Contrary to 
the phase separated state, the charge-density wave could 
survive the presence of long-range Coulomb interaction. 




FIG. 9. Real and imaginary parts of the spin-wave energy 
w q (U/t = 10 and x = 0.04). 




. 5 



FIG. 10. Spin-wave energy in the direction q = (q, q) 
for U/t = 10 and x = 0.04. Solid line: Re(w q ) (note 
the Goldstone mode at q — Q ~ 2.74), dashed line: 
Imfo/ q), d ot-dashed line: excitation spectrum at half-filling 
[Eq. Q3.4q )]. The dotted lines indicate the boundaries of the 
particle-hole excitation spectrum. 
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FIG. 11. Inverse propagators (pqP^) u=0 (solid line) and 
(qgqs)Z=o (dashed line) in the direction q = (q, q). 



When x ^> t 2 /U 2 (i.e. x/p 2 <C 1), it is tempt ing t o 
neglect altogether the x dependent terms in Eq. (3.44). 
This leads to 



u 



J(l-P 2 ) 
1-x 



X 



2 "1 2 



y q -2 



1/2 



(3.47) 



which is the result obtained by Arrigoni and Strinati 
within an RPA analysis in the regime t 2 /U 2 <C x <§C l.t£l 
Although Eq. ( |3.47 ) reproduces the main features of the 
spin- wave-excitation spectrum (see Fig. 1 in Ref. fl0| ), 
it incorrectly predicts a vanishing of w q along the line 
1x = %■ Even if there is no contradiction with Gold- 



stone theorem, the approximation leading to Eq. (3.47) 
is too crude to yield a correct description of the Gold- 
stone modes at ±Q. 



D. Spin- wave modes in the ferromagnetic phase 

When x > x c , the ground-state is ferromagneti c wit h 
p = 1, i.e. Q = (27r,27r) = (0,0). From Eqs. ( pp| ), 
fl02j ) and we obtain the spin-wave excitations 



(1 - x)^ = 2t(x - x c ){2 - T q ) + xJ(2 - Yq) 2 . (3.48) 

Since F q= o = 2, the mode at q = is gapless, as expected 
for a ferromagnetic spin arrangement. For q — > 0, 



J, 



off 2 

2~ q > 



JeS = 2t- 



1-X 



(3.49) 



The instability of the ferromagnetic phase against the 
formation of a spiral phase is signaled by a softening of 
the spin-wave modes, since the effective ferromagnetic 
exchange constant J c ff va nishes whe n x — > x c . Note 
that for x = x c , Eqs. (3.44) and (3. 48) coincide and give 



xJ{2 - Y q ) 2 /(1 - x). 
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IV. CONCLUSION 

In this paper, we have analyzed the semiclassical limit 
of the 2D Hubbard model in the framework of the spin- 
particle-hole coherent-state path integral. The main 
characteristic of the latter lies in the clear distinction 
between charge (7) and spin (SI) degrees of freedom. 
The semiclassical analysis consists in expanding around 
a broken-symmetry ground-state by making a saddle- 
point approximation on the spin variables, while the 
fermionic degrees of freedom are integrated out within 
a systematic t/U expansion. This should be contrasted 
with the Hartree-Fock/RPA theory whose validity in the 
strong-coupling limit is not obvious even at the semi- 
classical level. Nevertheless, we have shown that in the 
low-doping limit, the Hartree-Fock/RPA theory does re- 
cover the correct t/U expansion. We have also justi- 
fied the strong-coupling functional integral approach of 
Refs. where the electron-electron interaction is consid- 
ered within a large-C/ Hartree-Fock approximation while 
the SU(2) spin-rotation invariance is maintained by intro- 
ducing a fluctuating spin-quantization axis. We expect 
however these approaches to break down with increasing 
doping since they do not capture all processes of order 
t/U. In the spin-particle-hole coherent-state path inte- 
gral, a correct treatment of all processes of order t/U 
follows from the consideratiop-.of the quartic term of the 
action S [7*, 7; SI] [Eq. (§J)].El 



The spin-particle-hole coherent-state path integral pro- 
vides a convenient framework for studying spin-wave 
modes around the semiclassical ground-state. In particu- 
lar, it allows a systematic spin conserving analysis which 
ensures a proper description of the Goldstone modes of 
the system. From the mean-field fermionic self-energy, 
one can generate all vertex corrections to be taken into 
account in the propagator of the spin-wave modes. A 
mere calculation to order 0{t/U) of the effective action 
S[p,q] of spin fluctuations would give a contradiction 
with Goldstone theorem. Given this observation, it is not 
surprising that a t/U expansion of the RPA equations for 
the collective modes in the spiral phase meets with dif- 
ficultiep-jegarding a correct description of the Goldstone 
modes .E9 

Our main result is an analytical expression of the spin- 
wave excitations of the spiral phase over the entire Bril- 
louin zone except in the vicinity of q = 0. Besides the 
gapless mode at q = 0, we find two Goldstone modes at 
±Q as expected for a spiral spin order. Our result agrees 
with previous numerical calculations and improves the 
analytical expression obtained in Ref. fo[ Although the 
spiral phase is unstable, as shown by a negative mean- 
field susceptibility and the presence of imaginary modes, 
we expect the short-wavelength fluctuation modes (with 
real energies) to survive in the actual ground-state of the 
system. 



APPENDIX A: FREE ENERGY IN THE MEAN-FIELD STATE 



In order to calculate the free energy F [Eq. (2.17)], it is convenient to introduce a "Fermi wave vector" K^(a) = 
Kp (a) (cos a, sin a) defined by £k f («) = 0. From Eq. ( p. 16 ), we deduce 



K 2 F {a) = 2 



Atp- jU 



cos* a 1 sin a 



+ 



The free energy can be written as F = N 1 ^Z k (l — $(ek))ek where 



N 4-^ 



H- J(l- P 2 ), 

K F {cc) dR _ 



N 

K 

2 "da 
2n Jo 
(tip - m) 2 

47T 



2tt 



^m x m y . 



-K 



-fj, + Up ■ 



K I cos a sin a 



(Al) 



(A2) 



(A3) 



Eq. (2.17) follows from (A2) and (A3). 



APPENDIX B: CORRECTION F' TO THE FREE ENERGY 

The contribution F' to the free energy (Fig. ||) is given by 
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F ' = J^32 E E T U,U iu >'-> tuJ > iJ ) \ G f Ml^T (k, iw) (2*p cos k v f\Gf {iJ)\ 2 Q x (k', «/) (j2t sin | sin 



k,k' 



Performing the sums over cj and a;' and retaining the leading order term in t/U, we find (for T — > 0) 
F ' = ^E^kX 1 -e(e k 0)(2tpcosfc,) 2 (2tsin|sinfc^ 2 ! 



Using 



k,k' 



l^^( ek )4tV (C ° Sfcl/) 



N 



^E< 



(4ip) 2 



(e k =o + M) 2 



(Bl) 



(B2) 



(B3) 



we obtain F' = xJ(l — p 2 ) in the limit x -C 1. 
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